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The Bohr-Sommerfeld quantization rule and the Gamow formula for the width of qua- 
sistationary level are generalized by taking into account the relativistic effects, spin 
and Lorentz structure of interaction potentials. The relativistic quasi-classical theory 
of ionization of the Coulomb system {Vq = —^/f) by radial-constant long-range scalar 
(Si,r. = (1 — X){ar + Vb)) and vector (Vj.r. = A((7r + Vb)) fields is constructed. In the 
limiting cases the approximated analytical expressions for the position Er and width 
r of below-barrier resonances are obtained. The strong dependence of the width F of 
below-barrier resonances on both the bound level energy and the mixing constant A is 
detected. The simple analytical formulae for asymptotic coefficients of the Dirac radial 
wave functions at zero and infinity are also obtained. 
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1. Introduction 

Wide range of problems from various fields of physics (elementary particle physics, 
nuclear physics, physics of atomic collisions, etc.) is related to representations of 
formation and decay of non-stable (quasistationary) states of quantum systems.!^ 
Properties of such states are of interest for investigation of ionization of atoms, 
ions and semiconductors under the influence of constant and homogeneous electric 
and magnetic fields,'^ for description of cluster decays of atomic nucleP and effects 
of the spontaneous creation of positronSH, in consideration of a vacuum shell of 
supercritical atom,l2HIl and from the point of view of studying the Dirac equation 
in the strong external fields as well. 

The relativistic theory of decaying (quasistationary) states is elaborated quite 
well for the cases when components of an interaction potential of a fermion with 
external fields belong to the vector type, i.e. are Lorentz-vector components. 



27, 2010 13:56 WSPC/INSTRUCTION FILE Lazur 



2 V. Yu. Lazur, O. K. Reity & V. V. Rubish 

the same time the interior logic of development of the decaying states theory 
obviously dictates the statement of qualitatively new problems having an origin in 
nuclear physics and elementary particle physics. So, for example, quark structure 
of nucleons and multi-quark states in a nucleus compel to look at the nature of the 
inter-nuclear forces in a new way as well. Under the stimulating influence of QCD 
the problem of manifestation of quark-gluon degrees of freedom in atomic nuclei 
and nuclear processes has gained new development in the last decade and now, 
undoubtedly, constitutes the main prospect of fundamental researches in this field 
of nuclear physics. Without going into the further discussion of these important 
problems, note only that there are already enough reviews, books and papers where 
the specified tendency is sufficiently refiected (see, for example Ref.|Hl and references 
therein) . 

As is known, the problem of penetration through potential barriers underlies 
a theoretical view of the phenomena of cluster decays of atomic nuclei and non- 
stable resonance states of strongly interacting elementary particles. The subsequent 
theoretical description of such phenomena should be based on the relativistic wave 
equations by taking into account that besides electromagnetic forces the interactions 
between elementary particles can be also realized by the forces which are not de- 
pendent on the electrical charge. Take for example nuclear interactions between the 
nucleons, caused by interaction of nucleons with the meson field. One more class 
of (this time long-range) forces related to the considered problem are the forces 
arising between nucleons during exchange of electron-neutrino pairs (so-called (3- 
forces). For the first time such forces have been introduced by I. TamrrP in 1934, 
and more recently the vector and pseudovector variants of these interactions were 
considered in detail by authors of the well-known monography.^^ From the point of 
view of new problems of the strong interaction theory it is interesting to explore the 
more general case when a spin- 1/2 particle interacts with the scalar and vector fields 
simultaneously. As is known now there is a reason to think that such interactions 
exist between quarks in hadrons. 

The main difficulties of the theory of quasistationary states (for the applications 
specified above) are caused by the fact that in many cases the interactions cannot 
be described within the standard methods that use expansions in small energy 
parameter. In the problems related to the description of quasistationary states of the 
relativistic composite systems, additional difficulties arise when solving the Dirac 
equation with unseparable variables. In the modern theory of decaying states these 
difficulties are overcome by means of the relativistic version of the imaginary time 
method developed in papersBii£l This allows to calculate the tunneling probability 
of relativistic particles through potential barriers including those that do not possess 
spherical symmetry. 

Though this method has heuristic force and physical clarity, nevertheless it can- 
not be considered to be strictly mathematically justified, despite some attempts 
that were made in this direct ion EE] As is known, the accounting of the Coulomb 
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interaction between the outgoing electron and atomic core within the imaginary 
time method encounters considerable difficulties and, for example, in the theory of 
niultiphoton ionization of atoms^ is not fully carried out up to present time. 

Fortunately, many interesting questions of the relativistic theory of quasistation- 
ary states can be elucidated on an example of the simple models permitting exact 
or asymptotically exact solution of the Dirac equation. From the variety of such 



problems here we consider the hybrid version ( 13 ) of spherical model of the Stark 
effect (SMSE). Inclusion in the standard SMSE of "new" interactions, related to 
the scalar field, opens new possibilities for its applications in the relativistic nuclear 
physics and QCD. In a more general context we consider the nonstandard modeling 
problem of study of simultaneous influence of the radial-constant scalar and vector 
fields on a system of Coulomb levels. 



The quasistationary solutions of the Dirac equation in the composite field ( 13 ) 
atl/2<A^l and the corresponding complex spectrum of energy are generated 
by the radiation requirement meaning that at infinity the solutions F(r), G{r) 
are the divergent waves (see Ref. J2). If E is real, such solution does not exist, 
however, there are infinitely many complex quasistationary levels E = — iT /2 
with exponentially small imaginary part F. 



For the field ( 13 1 consisting from the mixture of Coulomb and radial-constant 
vector and scalar fields, the Dirac equation is separable in spherical coordinates. 
Apparently, this circumstance should essentially facilitate calculation of the position 
Er and width F of resonance. Let us note, however, that the usual approach (the 
numerical solution of the Dirac system) encounters the known difficulties related to 
exponential increasing of the Gamow wave function (at r oo) of a quasistationary 
state. In view of complexity of this problem we shall solve it in the quasi-classical 
approximation that gives the useful analytical expressions for the position E^ of 
resonance and its width F. Except the reasons of convenience the WKB method, 
or quasi-classical approximation, possesses a number of the principle advantages 
as compared to other methods. As is known, unlike the perturbation theory this 
approach is not related to the smallness of interaction and, consequently, has a wider 
range of applicability, allowing to investigate qualitative regularities in behavior 
and properties of quantum mechanical systems. Other important advantage of this 
method is in its applicability to the cases of both electromagnetic and scalar external 
fields. Further we shall consider the version of WKB method offered in Refs. IT2l 
1131 that can be used in the case of both the discrete spectrum and quasistationary 
states (resonances). 

This paper is organized as follows. In Sec. 2 we generalize the BohrSommerfeld 
quantization rule to the relativistic case where a spin- 1/2 particle interacts with 
scalar and electrostatic external fields simultaneously. In the following section by 
means of WKB method we solve the problem of finding asymptotic coefficients of the 
wave function in zero and at infinity with appropriate accuracy. Sec. 4 is devoted to 
the construction of the quasi-classical theory of ionization of Coulomb system {Vc = 
—£,/r) by radial-constant scalar (S'j.r. = (1 ^ A)((T7- + Vq)) and vector (V; = X{ar + 
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Vq)) fields by taking into account the relativistic effects and fermion spin. In the 
limiting cases ajim} <^ 1 and aj/E^ <C 1 the approximate analytical expressions 
for width of below-barrier resonances F which show the strong dependence of F on 
both the energy of bound level Er and mixing parameter A are obtained. 



2. Position of quasistationary states 

Separating the angular variables in the Dirac equation with the spherically sym- 
metric vector V{r) and scalar S{r) interaction potentials, we obtain the system of 
first-order ordinary differential equations for the radial wave functions F and G 
(c=l): 



dF k 

h^ + -F- [{E - V{r)) + (m + S{r))] G = 0, 
dr t 

h^--G+ [{E - V{r)) - (m + S{r))] F = 0. 
dr r 



(1) 



Here h is the Planck constant, k — hk, k — ^{j + 1/2) is the motion integral of the 
Dirac particle in a central field, jh is the total angular moment; the definition and 
normalization of the function F and G are the same as in the papersli2IH| 

For the description of the phenomena related to the formation and decay of 
quasistationary states, we consider the class of potentials V{r) and S{r) for which 
the effective potential (see ( [TT| )) of the squared Dirac equation possesses a barrier 
(of the type shown in Fig. [l|^ 

Let us give the algorithm of construction of quasi-classical solutions of the Dirac 
system As usually, we seek a solution of this system in the form of product of 



U(r,E) 




Fig. 1. Barrier-type effective potential U{r,E); ro, ri and r2 are roots of the equation p^(r) = 0. 
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rapidly oscillating phase factor and slowly oscillating amplitude: 

^^^^^ (gIo) "*^^^|/ j^^^')^*"'! "^W' 

oo 

y{r) ^ h-^y_i{r) + Voir) + hyi{r) + . . . , (^(r) = ^ /i".^^") (r), 

where /i — is a small parameter, (p{r) and ip^^^ (r) are the two-component quanti- 
ties (the upper and lower components of lyS*^"^ (r) correspond to the radial functions 
F and G respectively). Expansion of y{r) and (f{r) in h leads to a chain of the ma- 
trix differential equations for y(„)(r) and ip^"'\r) which are solved consequentially 
by means of known technic of the left and right eigenvectors of the homogeneous 
system. Let us give the final formulae for the wave function of a quasistationary 
state. These formulae have a different view in three regions: i) a potential well 
rQ < r < ri] ii) below-barrier region ri < r < r2', iii) classically allowed region 
r > r2 with quasi-discrete energy spectrum. Here rg, ri and r2 are the turning 
points, in which the radial quasi-classical momentum 



p{r) = ^{E - V{r)Y - (to + S{r)Y - {k/rY 
becomes zero (see Fig. [T]). In the classically allowed region tq < r < ri we have: 

E -V + m + S"'^'^ 



G(r) = C±sgnfc ~ ^^'^^"^ ~ ^ ^^^ cos 62. (2) 
Here we use the following new notations 

r r 
^ / N f f kw\ , TT ^ , . f f kw\ , TT 

0iM=y (^'+— j'^^+I' 02(r)=y jdr+-, (3) 

ri ri 

- 1 / V'-S' _1\ . _ 1 / V' + S' 1\ 
'^~2\m + S + E-V r)' ' 2\m + S - E + V ^ r ) ' ^' 

If a width of quasistationary level F <C 1 (which is justified by the answer) the 
condition of normalization by a single particle localized in the region < r < ri 
defines normalization constants C^: 

-1/2 

{E^^G^)dr=l^ \C-\ = {l^^^dr) .(|)''^ (5) 

where T is the period of radial oscillations of a classical rclativistic particle localized 
in the region I (rg < r < ri). 
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In the below-barrier region (ri < r < at fc > the solution, corresponding 
to the decreasing exponent, is of the form 



X 



G 



VqQ 



exp 



{m + S)V' + {E-V)S' 



2qQ 



dr 



m+S-E+V 



,(6) 



and for states with fc < 



X 



C2 

VqQ 



exp ■ 



(m + S)V' + {E- V)S' 




(7) 



Here q = \p{r)\, Q = q + |fc|r~^. 

At last, at r > r2 the divergent wave corresponds to the quasistationary state; 
at fc > the quasi-classical formulae for F and G have the form 



VpP 



exp ■ 



ip - 



{m + S)V' + {E-V)S'' 
2^ 



dr 



iP 

m+S-E+V 



, (8) 



where P —p + i\k\r ^. For states with fc < the radial wave function x is given by 



X = 



VpP 



exp ■ 



ip - 



{m + S)V' + {E -V)S'' 
2^ 



dr 



m+S+E-V 
iP 



(9) 



The quasi-classical representations ([2])- ([9]) constructed are invalid in small neigh- 
bourhoods of turning points rj {j = 0, 1, 2). To bypass these points and match the 
solutions one can use the Zwaan methocP^that allows to establish relations between 
normalization constants: 



G^ 



-tGi = T 



E-V{ri)+m + Siri) 



-I ±i 



\k\ rr 



X exp ■ 



q± 



(m + S)V' + (E- V) S' 

2^g 



dr 



(10) 



The quasi-classical formulae obtained above allow to solve a wide range of prob- 
lems in the theory of quasistationary states. So, in the case of exponential smallness 
of the barrier penetrability in the effective potential (EP) 



F — 

U{r,Er)^^V + S+— -(-— , 

m Zm Zmr^^ 



[Ei - m^)/2TO (11) 



the real part of level energy i?,. = En^k is defined by the quantization condition 



P 



; w 



pr 



dr = 



rir 



= 0,1,2, 



(12) 
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where rir is the radial quantum number. 

Now we explore one concrete example of the vector and scalar potentials 

V(r) = -^+\v(r), S(r) = (I- \)v(r), v(r) = <7r + Vo, (13) 
r 

in whose mixture at 1/2 < A ^ 1 and any a ^ quasistationary states of the 
Fermi particles exist; here Vq is a real constant, ^ is the Coulomb coupling constant, 
and A is the coefficient of mixing between the vector and scalar parts of the long- 
range potential v(r). The potential v(r) includes contributions of Lorentz-vector 
^.r.(^) — ^f^if) and Lorentz-scalar 5; ,..(^) = (1 ^ X)v{r) components, of which (at 
1/2<A^1) the first one dominates in all range of r, < r < oo. 



Even in the case of the relatively simple model of an interaction ( 13 1 the effective 



potential U {r, E^) (11) is a complicated function of the initial potentials V{r) and 
S{r), level energy Er, total moment j, and has an essentially different shape at 
^ A < 1/2, 1/2 < A ^ 1 and A = 1/2. The appearance of the quasistationary 



states in the potential model (131 at 1/2 < A ^ 1 and ct 7^ is very natural. Under 
specified conditions at small distances (r < 1) EP U (r, Er) corresponds to repulsion 
proceeding to the quadratic attraction at r ^ 1. It indicates at the definitive role of 
the long-range vector forces (electrostatic ones for example) in a barrier occurrence 
in EP U {r,Er). This ensures a possibility of level decay by particle penetration 
through the potential barrier in EP U{r,Er). Thus, at 1/2 < A ^ 1 the bound 
state is transformed into the quasistationary one. At r 00 its wave function has 
an asymptotic behavior of of a divergent wave type ([8|, At the same time, at 
any a ^ the energy E{a) becomes complex (E — Er — iT/2), and its imaginary 
part is directly connected to the probability of level ionization by the external field 



(13) 



Interest to the potential model ([T]), ( [T3| at 0^ A < 1/2 is basically caused 
by its use for the description of mass spectra of mesons and baryons and is not 
exhausted even till now (see, for example Refs. [T2l [TSl and references therein). It 
successfully combines the effective Coulomb interaction at small distances r with 
the linear increasing at large r that leads to a confinement of color quarks and 
corresponds to the string "stretched" between quarks. In our recent papei^ the 
detailed calculations of spectra of heavy-light mesons were carried out by means 



of the potentials (13). In this case the model parameters (13) are determined by a 
standardly: a is the string tension, ^ is the Coulomb coefficient, Vq is eigenenergy 
of a static source, and parameter A (0 ^ A < 1/2) gives the relative contribution of 
the vector VJ.r. (?■) and scalar S'z.r. (f) potentials. 

The term (5^ — V'^)/2m in EP ( 11 ) dominates at large distances and at ^ A < 



1/2 leads to the effective repulsion regardless of the sign of the parameter a. This 
repulsion is purely a relativistic effect and related to the fact that the interaction 
of fcrmion with the scalar external field S{r) is added to the scalar quantity m 
(particle mass), whereas the vector potential V{r) is introduced into the free Dirac 
equation by the minimum way as the temporal component of the 4-potential A^. 
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The presence of the relativistic terms (5^ — y^)/2m in EP U (r, Er) leads to the fact 
that U {r,Er) (1 — 2A)((Tr)^/2m) mcreases quadratically at large distances when 
^ A < 1/2. This ensures a confinement of color quarks and gluons (not observed 



in the free state). These properties of EP U{r,Er) (13) of the interaction model 
indicate that namely scalar fields, not vector ones, are essential for confinement. 
Considering the inter-quark potential, it is worthwhile to investigate the qual- 



itative regularities of the quasi-classical spectrum of the interaction model ( 13 ) in 
the phase of outgoing color, 1/2 < A ^ 1. One can even say that without such study 
it is impossible to completely understand the physics of confinement. 

Of course, the considered physical objects (quarks) cannot be free particles. 
However, it is not necessary to imagine that the string between quarks can become 
infinitely long. Outside of applicability of perturbation theory, at large distances the 
strong coupling inevitably leads to the creation of color charged particles and the 
problem becomes multiparticle in essence. Therefore, the problem about scattering 
of quarks, forming such colorless objects as hadrons, is usually considered. The same 
is valid for QCD problems related to study of quasistationary states of potential 
type. Here the problems of influence of the most various physical factors such as 
external fields with the mixed Lorentz-structure of interaction potentials, adiabati- 
cally slow variation of interaction parameters, connections with other (exterior and 
interior) degrees of freedom, translation and spin symmetry etc. on properties of 
decaying states arise quite naturally. 

Basically, one can not exclude that the results obtained when solving similar 
problems may appear useful for the description of tunnel effects in physics of con- 
densed matter (for example, in two-zonal semiconductors). So, for example, the 
electron in its motion in a solid behaves as a particle with "effective" mass JTie//, 
which differs considerably from the true mass m. Moreover, We// can depend on 
distance r and motion features, because the difference of masses {meff — m) is ac- 
tually caused by interaction of an electron with objects surrounding it. In the same 
sense all quarks masses are "effective" because they are determined in processes in 
which a quark interacts with other quarks. 

In the following two sections we shall set out the applications of elaborated quasi- 
classical formalism to calculation of shifts and widths of quasistationary states of 
the potential model (13 1 with outgoing color charges (1/2 < A ^ 1). Special case 
A = 1/2, separating a phase of confinement of color charges and a phase of outgoing 
color (i.e. deconfinement), we will consider in our next papers. 

Following the described in Ref. [TH (see Sec. 4 therein) calculation scheme, we 
represent the quasi-classical momentum in a convenient form: 

Pir) . I.IV2A3TM1 = HV2A3I[("--)(^-)(--'-)(---^)]^^^ 



Here cr plays a role of intensity of the radial-constant long-range field, and the 
turning points a, 6, c and d are solutions of the equation r'^ + fr'^+gr'^ + hr + l = 
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with coefficients 



/ = 



2??i 
(1-2A)ct' 



9 = 



(1-2A) (72 ' 



h = - 



2Eri 
(l-2A)cr2 



/ = 



(1 



r_ 

- 2A) a2 ^ 



where 7 — -^/fc^ — ^2 g^j^j _ (1 — A) m + A-E^- The characteristic parameters 
Er = Er — AVo and m = m + (1 — A)Vo introduced here mean the "shifted" energy 
and the "shifted" mass, respectively. The specified equation has four real roots 
(d < c < 6 < a), determined by the expressions 



where 



a, b 



4 



4 



2g 
3 



1 



(5±Ah 



c, d 



1/3 



2 



4g 
3 



D = -f +4:fg-8h, u = -3fh + 12l. 



1/3 



1/3- 



1/3 



/ 
4 

1/2 



1 



:sta_), 



(14) 



4H' 



Z = v+ y/-4:U^ + v'^, 



= 2g^ - 9fgh + 27h^ + 27 f I - 72gl, 



the upper signs in ( 14 ) correspond to the turning points a and c, and the lower ones 



correspond to the points b and d. 

In the quantization condition ( |l2[ ) we integrate over the range of r where E — 
U{r,Er) > 0. For the potentials ( |13[ ) considered by us at 1/2 < A ^ 1 this means 
that vq — c and ri = b where c < r < b and r > a are classically allowed regions, 
and b < r < a is below-barrier region, in which p^(r) < 0; at r > a the particle 
goes to infinity. If the maximum of the potential function U{r,Er) is at the point 
Tmax « 77i[(2A-l)cr]~^ for tj > (jmax ~ [£^r^/('7i kl )] for 0" < 0), and minimum 
is at Tmin ~ l"^ I ErS, then the energy spectrum of below-barrier resonances is in the 
range Umin < Er < Umax ■ The considered situation is schematically shown in Fig. [l] 
where black dots indicate the position of turning points. 

Using the technique of evaluation of quantization integrals from Sec. 4 of Ref.jT^J 



we can represent the phase integrals of the equation ( 12 ) in the form 

b b 
Ji= I p{r)dr V2A-l|cr| 



(r^ + fr"^ + gr + h + Ir^^) 



Ri 



dr. 



(15) 



J 2 = 



k 1 



-dr ^ — 



where 



p{r)r 2V2A - 1 \a\ 



A+ = 



dr 



2cr(l - 2A) 



dr 



{r-X+)Ri J (r-A_)i? 



E + mT\ nE + to)2 - 4Cct(1 - 2A) 



(16) 
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In the general case for the arbitrary values of ^ and a the integrals Ji and J2 cannot 
be expressed in elementary functions. The carried out transformation of Ji and J2 
is convenient because after replacement of the integration variabl^^ 

c{b — d) — d{b — c) sin^ Lp 



b — d — {b — c) sin^ Lp 



(17) 



the integrals from right hand sides of the formulae ( 15 ) and ( 16 ) are transformed to 
the complete elliptic integrals of the first, second and third kindll51i3 As a result, 
we have obtained the transcendental equation 



2V2A- 1 
V(a-c) (b-d) 



k 

^ 2(2A-1 



[NiF{x) + N2E{x) + N3Tiii>,x) 

[{c-d) {N,n{D+,x)+NeTl{i?_,x)) 

+7V7^^(X)]}- (rv + ^^TT, (18) 



defining (in the quasi-classical approximation) the real part E,. ~ En,.k of complex 
energy of quasistationary states at Umin < E^ < Umax- Here the quantities P±, 
X, ^, Nj {j = 1, 2, 7) are obtained from respective expressions (A.l |-( A.6) for 
^1 ^±1 Xi (recently found by us in Ref. I12p by making the simultaneous 

replacements a^b, b^c, c~^d and d ^ a. 

For the determination of the Stark energy Er = En^k the equation ( 18 1 can 



be solved by numerical methods only. To construct the asymptotic behavior of 
resonances it is necessary to apply asymptotic methods of calculation of phase 
integrals just as in Sec. 4 of Ref. [T^] This imposes some restrictions in calculation 
of shifts of quasistationary levels and its widths for both small values of intensity 
cr of radial-constant (scalar-vector) long-range field and not too large ones. Namely 
the quantity m divides the range of Er = Er ~ AVq into the two domains, in which 
the spectrum of quasistationary states has a various asymptotic behavior. Consider 
some of the most typical situations connected with the relative values of the energy 
Er and level to. 

Case A^If Er < rh, Umm < Er < m (see Fig. [l] where U,mn = C/(r™i„, E'^), 
Tmin ~ 1^ I ErC] ^-nd Condition a /S^fh? <C 1 is satisfied, the pair of classical turning 
points a and d is rather distant from pair of points c and 6 (see. Ref. I12p : 



Er^±( 



1 - 



Er^±^ 



m — Er 



m + Er 



o' m- Er 

Hereafter, we use the notation 



cr(l-2A) m + Er' 



(19) 



(20) 



= d [ErkY - {m^Y, ^=-Jm2~£;2, 772 = Am + (1 - A)^;,. 
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In this case the derivation of asymptotic expansions of quantization integrals (15 1, 
( 16 1 in a sniaU parameter cr/^m^ is carried out just as in item A of Sec. 4 of Ref. [T^] 



and gives the former expression for the level energy 



En 



2Cm2 



TO 



K+7) 



mo 



-1/2 




2£,'^mEo 
Mo 



n'^ = Ur + {I + sgnfc)/2, 



and the quantities fiQ, r]io and 7/20 are obtained from /i, rji and 772 by replacing Er 
hyEo. 

The influence of weak radial-constant scalar and electric fields on the system 
of Coulomb levels has been analyzed before (see, for example Refs. [121 fTSjl . The 



analysis was carried out both on the basis of quasi-classical formulae (21) and 



by the numerical solution of the transcendental equation (18 1. In particular, the 
calculations have shown that the sign change in ct (cr —> — cr) leads to small changes 
of energy spectrum, if \a\ <C 0.2 GeV^. 

Case B. Let us now consider fh < Er and m < Er < Umax (see Fig. [l] where 
Umax = U{rmax,Er) and Tmax ~ ?7i[(2A - l)cr]"^). Here the ratio aj/E^ plays a 
role of the small parameter. In this case the quasistationary states in the composite 



field (13) exist only at the positive values of parameter cr. 



Under the requirements formulated above, from ( 14 ) we have obtained 



c, d '. 



-Er£, ± e 



(23) 



Er + m ^ ^ Er — m 

~cr(2A-l) Er + m ^ cr Er-m ' E^ 

It is seen that turning points c and b are rather distant from one another (a, b 3> 



c, \d\). This allows to evaluate the quantization integral (12) analytically. 



Further we shall give only the recipe of evaluation of the quantization integrals 



Ji_2 . Just as in the item B of Sec. 4 of Ref. [121 the integration range in ( 15 ) and ( 16 ) 



we divide into two domains c ^ r ^ r and r ^ r ^ 6 by introducing the dividing 
point r satisfying the requirement c ^ r ^ b. In the first domain c ^ r ^ r we 



calculate the integrals ( 15 ), ( 16 ) by expanding the quasimomcntum p(r) in a power 
series in the parameters r/a^l and r/b<^ 1. In the second domain r ^ r ^ b the 
expansion of p{r) we carry out in the small quantities c/r <^ 1 and \d\/r <C 1. 

When we add the asymptotic expansions of integrals over c ^ r ^ r and r ^ r ^ 
b the final result will not contain the quantity r. So, we have obtained the equation 
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If we expand the left hand side of (24 1 in ffi/Ej. <C 1 up to the terms of the third 
order, then we obtain the transcendental equation for the level energy which we 
solve by the method of consecutive iterations. Thus, we arrive at the expression for 
the energy (within 0{a^ / E"^)) 



2a{l - 2A) ( e log ^'^'i^ + 3e + KA + -kN 



4£;(o) ■ 



+\m'^{l - XA) 



1/2 



+ Ayo, (25) 



where 



£;(o) ^ ^(0) _ xVo, A^{2\- 1)-^''^ log [(l + A + \/2A^) (1 - A)-i] , 
C= (l-A)2A-A-2(7e(l-2A)/(i?("))^ B = (l-A)(l-Al)m-4(7^(l-2A)/i?(°\ 



1 sgnfc If f i \ ^ 

A^ = -.+ 2 + — + -(^7arccos 



(2A-l)-i/2log 



77i + V(2A-l)(i?2 



(26) 



and = i?'-"^ — AVq, is the zeroth approximation of energy, by the choice 
of which the quantity En^k depends very weakly and in most cases one can take 

^(0) wi;,(e = o). 



The expressions ( 24 ) , ( 25 ) differ from the corresponding formulae in the case 
of purely discrete spectrum^ only by replacements rj ^ rj and A — > A. So, the 
substitution rj ^ rj, A ^ A transforms the equations E^ 



WKB(as) 



obtained in the 



case of purely discrete spectrum, into the equations for a real part Er of energy ( 24 1 , 



(25 1 of below-barrier resonances. Accuracy of the calculation of E^. by the means of 



formulae (24), (25) is fully appropriate (see Sec. 4 of Ref . [T^ . and usually there is 
no point to make the result more precise for practical purposes. 



3. Asymptotic coefficients of a wave function 

Asymptotic coefficients Cf.g a-t the zero and Ap^c at the infinity are the character- 
istic parameters of a wave function. We shall give simple analytical approximations 
for these coefhcients which describe the results of numerical calculations quite well. 

First of all we consider the construction rules of asymptotic expansions of solu- 
tions of the Dirac equation at zero (r — 0) along with more standard expansions of 
solutions, when r — >■ oo. For the considered potentials (18) we have 



F,G = CKGr^ + ..., r^O, Cp/Ca = {k - ^)/^, (27) 
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and for wave functions of the discrete spectrum (0^ A < 1/2) the normahzation 

oo 

condition J (^F^ + G^) dr — lis satisfied. Values of F^{0), G^(0) (or, more precisely, 



of Cpq) define the probability to discover particles at small distances from one 
another and are of considerable physical interest especially in the case of systems, 
in which interactions of two various types (for example, the Coulomb interaction 
and long-range one v{r)) exist. 

In the classically forbidden range ^ r < tq the wave function of oscillating 
type ([2]) is changed by the solution decreasing exponentially with increasing r (see 
Fig. [l|T Matching the WKB-solutions of the Dirac equation on both sides of the 
turning point tq, for radial distribution functions F{r) and G{r) we obtain quasi- 
classical expressions in below-barrier range r < tq: 



F(r) = (-l)"- 



G{r) = sgnfc(-l) 



nfc(-l)""^ 



E-V + m + S 



q{r) 

E -V -m-S 
q{r) 



1/2 



exp 



1/2 



exp 



q ) dr 

q r 



q r 



, (28) 
• (29) 



All integrals in (281 and (29) are expressed through the quite complicated combi- 
nation of the elliptic integrals. But in the cases Ej. < rh and E^ > fh they can 
be calculated through elementary functions, using the relations a/^fn^ <C 1 and 
aj/E^ <C 1 to expand the integrands into power series. 

Let us first investigate the asymptotic behavior of the quasi-classical solutions 



(281, (29) at r for the lower levels {E < in) which are basically defined by the 
Coulomb potential (ajiw? ^ 1). Note that the larger the Coulomb parameter ^, 
the smaller is the essential potential v(r) at small distances. Before the evaluation 



of the asymptotic coefficients Gf,g by means of formulae (28), (29) it is necessary 



to expand the quasi-classical momentum p(r) in potential v(r). Then, using the 



technique of evaluation of the phase integrals from Sec. 2 and proceeding in ( 28 1 , 



(291 to the limit r — >■ 0, we obtain in zeroth approximation the expressions for the 



asymptotic coefficients at zero: 



|Cf| = 




iE^ + 



MO 



(30) 



Here Qq - 
formula 



[Eq fc)2 — (to 7)2, and the period of radial oscillations T is given by the 



(31) 



January 27, 2010 13:56 WSPC/INSTRUCTION FILE Lazur 



14 V. Yu. Lazur, O. K. Reity & V. V. Rubish 



When we derive the expression ( 30 ) we use the quasi-classical requirement ol the 



normalization Solving the Dirac equation ([ij at small distances (in range < 
r < b) one can neglect the term with the linear potential {a — 0). Having used 
the asymptotic behavior of the normalized radial functions F(r) and G{r) of the 
relativistic Coulomb problerrfl^ at r — )■ and the relations (27), we find the more 



exact (than pO| )) expression for Cp- 

{m + ^o)r(27 + < + 1) 



7+1/2 



(2mo) 
r(27 + l) 



Mo 



i' I 



1/2 



(32) 



Here = + (1 + sgnfc)/2. The formulae (30 1 and (32) differ one from another 
within an error between the Stirling formula 

n! = V27rexp{(n + 1/2) logn - n} (1 + 0(n"^)) at 



->-oo 

and the F-function. 

For states with Er > fh, when the requirement aj/E^ ^ 1 is satisfied, the 
Coulomb potential is essential only in the range of small distances and can be 
considered as a small perturbation in the basic range of particle localization (i.e. in 
classically allowed range b < r < a). This gives the possibility to exclude {a = 0) 
the linear part of the potential v(r) from the quasi-classical momentum p(r) when 
evaluating the integrals in exponents (28), ([29|. Then at r — the asymptotic 



behavior of radial wave functions F(r) and G(r) obtained in this way allows to 
determine the asymptotic coefficients: 




^{jm+\k\Er) 



exp 



^Er 



^Er 



El^m-^ 



,(33) 



where quantity 9 is defined in (19), and energy E^ is given by the formula ([18 
Characteristic feature of the considered case is the fact that in the integral ([5 
which defines a period of radial oscillations T, only the range of values of the inte- 
gration variable r, where the Coulomb potential can be considered a perturbation, 
is essential. By neglecting the Coulomb interaction, we arrive at the expression 

2 



-\\^E'^-m-^ + {l-\)^r]2 



(34) 



a{l - 2A) 

The asymptotic coefficients Ap^ Aq of radial wave functions at infinity are im- 
portant physical parameters of bound states as well. Along with the coefficients 
Cp, Cq at zero (27), the asymptotic coefficients ApQ are continually encountered 



in quantum mechanics,'^ atomic and nuclear physics)22Ei] in the inverse problem 
of quantum scattering theorjESEU etc. For the potentials (13) the quantities Ap^a 



are related to asymptotic behaviors of the normalized radial wave functions by the 
relations 

Vl - 2Acr 2 



F,G^ ApGT^ exp 



VI - 2A 



OO, 



(35) 
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where a > 0, Ap = -\J\ - 2A A 



G, 7 = 



2 
»?2 



A? 



, and the parameter 



2(l-2A)3/2a yr^2A' 
A has values in the range ^ A < 1/2. 

In the below-barrier range r > r\ = a far from the turning point r\ = a under 
the requirements ajtrn} <C 1 and E < fh, after the evaluation of integrals the 
quasi-classical solutions (|6|, ([t]) are of the form of decreasing exponents 




where = (1 — 2A)(r — c)(r — <£). The estimates show that by the satisfaction 
of the requirements a j^m} ^ 1 and E < fh there is quite long range of distances 
r which are much larger than size of the Coulomb hydrogen-like system (r ^ (r) , 
see (37) or (38) in Ref. IT^ and much smaller than the distance r w (i^^/r/icr)^/^ at 
which the Coulomb interaction becomes quantitatively comparable with the long- 
range interaction. In this range as the wave functions of zeroth approximation it 
is natural to take the unperturbed radial functions F and G of the relativistic 
Coulomb problem, and the potential v{r) can be considered as a small perturbation. 
Neglecting it, we arrive at the following quasi-classical expressions for F and G 







ym + 


Eo 


-\Jm 


— Eq 


\J fh + 


Eo 


-Jm 


-Eo 



A' 



m + Eo \ f £.m- k^ioY^^ ( ^Eq - 7/io \ ' f 2fi^r\ 



7/2 



(37) 



Equating (37) to the asymptotic (at r — cxi) representation of solutions of the Dirac 
equation in the Coulomb fielcP^ 




m + Eo \ . 



(38) 



' TO — i?0 

we obtain the explicit expression for the period of radial oscillations of the classical 
relativistic particle: 



T 



1 



~ \ 1/2 / 9\ / ~ \ T 



, £,Eo + 7/io , 



(39) 
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Here Ac is the asymptotic coefficient of the Dirac radial wave functions in the 
Coulomb potential: 

-,1/2 



\Ac 



2 ^ m2 p 



1 



(2/io) 



(40) 



Comparison of ( [37| and ( 38 ) shows that their exponential and power factors are the 
same, however, the asymptotic coefficients A^^^ and Ac differ within an error 
between the Stirling formula and F-function. 



The formula (37) is obtained by neglecting the linear part of the potential v{r). 
This approximation was argued above by the means of the circumstance that under 
the quasi-classical requirements (ct/^to^ <^ 1) there is a range of distances (r) <^ 
r <C r, in which the distortion of a wave function, caused by action of the linear 
part of the potential v{r), can still be neglected and there is the law of decreasing 



radial wave functions ( 37 ) that is characteristic for the relativistic Coulomb problem. 
Change of the law of dccrcascmcnt (37) of functions F{r) and G{r) to (35) at r ^ f 



appears because in EP U{r, E) we have taken into account the quadratic (in or) 
terms which increase with increasing r more rapidly than others and so play a role of 



the perturbation which destroys the asymptotic regime ( 37 ) . As a result of such an 



account, using the quasi-classical approximation (36) for the normalized radial wave 
functions F and G at large r, we obtain the following expression for the asymptotic 
coefficient at infinity 



Af = 2/ioAc(l - 2X)'<+^l^a'' 



A»0 



X exp 



2A^o + Vw 



4(1-2A)3/2ct 



(41) 



We now proceed to the other limiting case a-f/E'^ ^ 1 when the centrifugal 
potential '~f'^/2mr'^ does not play an essential role at large distances and can be 
omitted in the quasi-classical momentum p{r) = iq{r). Having expanded the quan- 
tity q{r) = \p(r) \ inpowers of the Coulomb potential and calculated the integrals in 
exponents in ([6|), ([7|) under the requirement a^/E"^ ^ 1, in the asymptotic domain 
r "> oo we arrive at formulae of type of (35) for F and G, in which 

(1-2A)1/4 /2(l-2A)a^ 

S,mri2 , iE j 771 

m, 



Af 



exp 



4(1 - 2A)3/2a 



2yr^2A(£;2 - m2) y/W^^ 



(42) 



and the period T is determined by the previous expression (34) 



4. Width of quasistationary states 

The quasi-classical approximation (or WKB method) is the most often used method 
of the approximate solution of the relativistic quantum mechanical problems and 
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leads to an obvious physical picture of percolation of a particle through the potential 
barrier in EP U{r,E). 

In the consideration above the quasistationary character of the Stark spectrum 
was ignored. Thereupon it is necessary to remind that for 1/2 < A ^ 1 and any 
value of a 7^ 0, U{r,Er) has the shape of a potential with a barrier, owing to 
what there are quasistationary states with complex energy E = Er — iT/2 instead 
of discrete levels. This is a result of the attraction of the term —V'^/2 which at 
large distances and 1/2 < A ^ 1 strongly suppresses the contribution of all other 



summands in ( 13 1 and transforms the finite domain of the fermion motion into the 
infinite one. Thereby one can assert that the infiuence of the long-range vector field 
V;.r.(r), dominating at large distances, is manifested not only by the means of a 
modification of the energy spectrum of the system, but also leads to a non-zero 
probability of its decay due to a fermion passing through the potential barrier in 
EP U{r,Er). 

The probability of tunnel transition of a particle from the bound state (with 
energy E^) into the continuum state is defined by the imaginary part (i.e. by the 
width F) of complex energy of quasistationary statesl^^ 

F = -2Im[G*(r)F(r)],^oo. 

Having calculated the fiux of the particles outgoing to infinity by means of quasi- 
classical formulae ([8])-(10), we find the following expression for the level width F 

r= iexp[-2!i], (43) 
T^2j'^,r. 0./(,-^).. (44) 



The obtained quasi-classical formula ( 43 ) is the relativistic generalization of the 



well-known Gamow formula for the width of a quasistationary level. The nontrivial 
moment of such a generalisation is the modification of expression for the period 
of oscillations T and the occurrence of the additional factor in the preexponent of 



expression ( 43 ) that depends on a sign of the quantum number k and is caused by 
the spin-orbit coupling in the mixture of the scalar S{r) and vector V{r) potentials. 

Thus, in the quasi-classical approximation the problem is reduced to evaluation 
of two characteristic phase integrals T and Vl. 

Having used again the notations of Sec. 4 of Ref. il2j we write the quantity q 



from ( 44 ) in a form convenient for our purposes 

r 

where the parameter A ranges values from 1/2 up to 1. By an appropriate (simi- 



lar to (15), (16)) transformation of integrands the integrals in right hand sides of 
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formulae (44 1 are reduced to the complete elliptic integrals. Without detailing the 



aforementioned calculations of T and il, we give only the final result 



dEr + £,- Xa \(f - 



WWia-c) [b-d) (2A~1) 



\(TV (c — dy 

m 



E{x) + {c~d) 



E~ Xa 2d 




Fix) 



n(^^,x) , (45) 



2V2A - 1 

V(a -c){b- d) 



NiF ix) + N2E (x) + iVsH (v, x) + 



2(2A-1)|ct| 



[ib ~ c) (TVsn {v+,x) + iVen [v^ , x)) + NrF (x)] 



(46) 



The quantities ly, y±, x-, ^1 ^±1 ^, {j ~ 1,2, ...,7) belonging to 



(451, (46) are defined in (A.1)-(A.6) and (18) 



The derived formulae (45), (46) are valid at large values of modules of phase 



integrals T, £7 and together with (43) solve the problem of calculation of width T 
of Stark below-barrier resonances at Umin < Er < Umax- However, these formulae 
are rather cumbersome and not too convenient for concrete calculations. With a 
purpose of deriving an analytical expression for the width of quasistationary level T 
the calculations for the cases Ej. < m and > m should be carried out separately. 

Case A. Let us begin with the simpler (in the sense of calculation) case of 
quasistationary levels with > fh [E^ < Umax) when the under requirements 
(T7 ^ E^, a > the classical turning points b and a are rather distant from 
the pair of points d and c. Asymptotic expansion of the barrier integral f2 can 
be constructed by means of the procedure which is very similar to the procedure 
applied to quantization integrals Ji and J2 in the case of purely discrete spectrum 
in the item A of Sec. 2. Omitting the details of the calculation, we give only the 
final formula for the width of the quasistationary level: 

r« lexp[-217(^„A)] , (47) 



n{Er,x) 



4 



2V2A - 1 I f7(2A - 1) 



2^A 



2Er^^/2X - 1 



Ei - m2 



(48) 



where 772 is defined above in (19). Corresponding expansions for the energy i?,. and 



period are given by previous formulae (24), (25) and (34). Analytical expressions 



(34), (47) together with the known dependencies of Er{cr), obtained by means of 



the formulae (24), (25) or by numerical solution of the transcendental equation 



(18), allow to calculate the widths of below-barrier resonances in "transitional" 



range of intensity a of the radial-constant (scalar-vector) long-range field. As can 
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be seen from the Fig.[2]the function fl{ErX) decreases monotonically with increasing 
the parameter A. Therefore, decreasing the relative weight (1 — A) of the Lorentz- 



scalar 5*;.^. in long-range part v(r) of the interaction (131 rapidly increases a 
probability of ionization of quasistationary level. 

When the level energy is close to the top of the barrier {Er Umax), the quasi- 



classical formula ( 47 1 for the level width becomes invalid. From the point of view of 
the effective potential it corresponds to closing of the turning points b and a owing 
to what the exponential smallness of level width F disappears. In this case in the 
below-barrier range b < r < a EP U{r,Er) becomes parabolic, and the quantum 
mechanical problem, similar to calculation of penetrability of a parabolic barrier, 
arises. The formulae necessary for considering this case are given, for example, in 
papersBilH^ 

Case 'B. At Er < rh, Umin < Er the asymptotic expansions of T and fl in pos- 
itive powers of small dimensionless parameter a j^rr? are constructed by the same 
technique, as in the item B of Sec. 2. In short, we shall derive only the asymp- 
totic approximation of the integral 17 which defines the barrier factor. Assuming 
that the potential barrier in EP U{r,Er) is rather wide, we divide the integration 
domain b ^ r ^ a into two segments by the point r* which satisfies the condition 
6 ^ r* ^ a (it is indeed possible because a — >■ cxi as u — >■ 0). In the first domain 
b ^ r ^ r* , the long-range potential v{r) can be considered as a small perturbation, 
and in the second one r* ^ r ^ a, on the contrary, the Coulomb field is much 
weaker than the long-range field and can be considered as a small perturbation. 
Expanding the quantity q{r) in small perturbation in each domain, we obtain some 



tabular integrals in ( 44 1 , the sum of which gives the value of the barrier integral fl 
to within terms Oia j^w?). By omitting the details, sense of which is clear, we give 
the asymptotic behavior of the imaginary part of energy of quasistationary state in 




Fig. 2. The function Q.{Er, ^) that defines the dependence of the exponential factor in the ion- 
ization probability (j47]l on the level energy Er, 0.44 GeV< Er < 2 GeV. 
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a weak-coupling regime: 



1 

f 



X exp 



|o-|772o6'o 

mi 



A»0 



m ^ 
1 



^0 k 



1/2 



9 



^ ^0 C + Mo 7 , 
2A ^ cr I arccos 



J?! 



sgn a 



.(49) 



ct(2A-1) |o-|V2A-1 V2A-1 / \ t]2 

This result is valid for both positive and negative values of parameter a, and in the 



potential well c < r < b one can again use the formula (31 1 or (39) for the period of 



oscillations of a classical relativistic particle with energy Er- If one uses more exact 



expression for the period (391, derived by matching WKB solution (37) with the 



asymptotic behavior of the relativistic Coulomb wave function at r — > oo, then one 



can represent the width of quasistationary levels (|49|) in the form 



r = 2^0 \Ac 
2sgna- 



(1 - A)??20P 

2A - 1 



exp ■ 



AC 



^(^o,A) _ 2Xfiap 
\a\ 2A-1 



VlQ _ 

mo 



-sgncr 



(50) 



where Aq is the asymptotic coefficient ( 40 ) of the normalized wave functio n in 



the Coulomb potential, and notations Eq^ /ig, 7710 and 7720 wer e in troduced in (21), 



(22). The functions <^{Eq,\) and p{Eo,X) from the exponent (50) are given by the 
formulae 



a>(i?o,A) = (2A-l)- 



{ \/2X - 1 



mo 

1 

mo^ 



-sgncr 



Vlo^J■o sgn a 



(51) 



2^m2 



Mo 



2^'^mEo 
Mo 



- k 



mo 



Comparison of the results of calculations of level energies based on the "exact" 



quasi-classical formulae ( 43 ) , ( 45 ) , and ( 46 ) with results of numerical calculations 
shows that relative error of (50) and (51 ) does not exceed 2% at \a\ ^ 10~® 



10- 



GeV^. So the formula (50) is convenient when used for quick estimates of F. In 
Fig. [3] the dependence of the function ^{Eq, A) on level energy_Eo are represented 
for several values of mixing parameter A. As is shown in Fig. (|3a, a < 0), $(i?0i A) 
increases with sinking of the level Eq and decreases when the mixing coefficient 
A (1/2 < A ^ 1) increases. The last fact can be easily explained: by increasing 
the relative weight A (from 1/2 up to 1) of the Lorentz-vector V/ ,..('') in the long- 
range part of [/(r, i?o) the attraction increases at large distances and the effective 
width of a barrier decreases. It results in the fact that at \a\ ~ 10~^ — 10~^ GeV^ 
the increasing of the parameter A (1/2 < A ^ 1) leads to the increasing of the 
ionization probability F and, on the contrary, decreases Eq or, in other words, the 
sinking of the bound level decreases F (at the fixed values A and a). In the case of the 
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positive values of a the function ^{Eq,X) decreases monotonically with increasing 
the parameter A (Fig. |3|, therefore by decreasing the relative weight (1 — A) of the 
Lorentz-scalar Si,r. (r) in the long-range part v{r) of the interaction ( 13 1 one rapidly 
increases the probability of ionization of quasistationary level (at the same value of 
a). 

The formula (50 1 is asymptotically exact at the limit — > 0. At small a ^ 



only the exponential factor, depending extremely rapidly on intensity of a long- 



range field and varying within many orders of value, plays an essential role in (50 1 



However, the applicability field of the formula (50) is rather narrow: \a\ < 10 



0,05 



0,00 



35 



34 



10 






X= 0.6 ^^^-^"^"^ 




X=0.1 




X=0.8 




X= 0.9 




J 1 h- \ 



0,94 



0,96 



1,00 



Fig. 3. The function ^{Eq, A) that defines the dependence of the exponential factor in the ion- 
ization probabihty (|50[l on the level energy Eq (in GeV): a) for a < 0, h) for cr > 0. 
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GeV2. 



As can be seen from (50), the width of quasistationary level is proportional to 
the squared asymptotic coefficient Ac- It is no wonder: at a/S^in'^ <C 1 and Er < m 
the ionization goes from the "tail" of the Coulomb wave function and the barrier is 
wide. 

In the case when in addition to the Coulomb field there is only the radial- 



constant electric field (A — 1), from the formula ( 50 ) at Vo = 0, (t < and ^ = Za 



{Z is the nuclear charge, a ~ 1/137 is the fine structure constant) we have obtained 
{h = m = c = 1): 



2^lo \Ac 



24 



exp 



2Za arccos Eq — 2pd 1 — Eq 



,(52) 



where Eq is the energy of a bound state in the absence {a = 0) of an external 
long-range field, and 



1 



P = 



2Za 



^ . 3Z^a^ ,2, , 
E,[^^-e]-k 



The function ^(Eq) in the exponent is given by the expression 



$(i?o) = arccos So - Eo \Jl - Eg, 
and possesses the obvious property $(— i?o) = tt — $(-Eo)- 



(53) 



Consider some of the limiting cases of the derived expression ( 52 ) : 



(1) Let us begin with ionization of s-level bound by short-range {Z = 0) forces 
under the influence of the radial-constant electric field a < 0. In this case from 



(52 1 we obtain the expression 



r oc exp 



(54) 



which coincides with the result of papers!^" for the Stark ionization of s-level 
bound by short-range potential or (5-potential (which is a good approximation 
in the case of ionization of the single charged negative ions, such as H~ , A^a" 
etc.) within exponential accuracy. 
(2) In the presence of the Coulomb fleld it is worthwhile to consider the various 



limiting cases for the quantities appearing in the formula (52): 



l-E, 



.1/2 



1 



arccos i?o = < 



l-E^n 



:Ei 



7r-(l-i?2 



2\l/2 



l-E, 



2\3/2 



Eo ^ 1, 

Eo^O, 
Eo ^ -1, 



(55) 
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r 2^/^ 



2 



20 



(1 - Eo) + 



25/2 



a+Eof^' 



, Eo^l, 
Eo^O, 
Eo ^ -1. 



(56) 



In the nonrelativistic limit {Eq — > 1, a — >■ 0) the formula (52) transforms into 
the known expressiorPSl for width of below-barrier resonances in spherical model of 
the Stark effect for hydrogen atom: 



r "^f-tji.r. I^n.r. 



2fi. 



2 

n.r. 



exp 



2^: 



3 

n.r. 



3|a| 



(57) 



where fJ-n.r. =^ Z/n, n is the principle quantum number, A^.r. is the asymptotic coef- 
ficient (at infinity) of the Coulomb wave function of the free (cr = 0) nonrelativistic 
hydrogen-like atom. 



The correction of the order of a in the exponent of ( 52 ) slightly increases the 



probability of ionization when compared with the corresponding nonrelativistic for- 



mula (57 1. The factor ^{Eq) in the exponent of (52) increases monotonically with 
sinking level (it is equal to it/2 and tt at E'o = and Eq = —1, respectively) which 
leads to rapid decreasing the probability of ionization. 

At Eq —1, that is for the level which has sunk to the edge of the nega- 



tive energy continuum, the leading (exponential) factor in ( 52 ) becomes equal to 
exp(— 7r/|fT|) and coincides with the corresponding factor in the Schwinger forniulaP^ 
(obtained within the quantum field theory) for the probability of creation of the 
electron-positron pairs from vacuum in the constant electric field. 

By means of the formulae obtained above the spectrum of quasistationary lev- 
els is described for the accepted hybrid version of SMSE. Such model qualitatively 
reproduces the following characteristic features of quasistationary states in an mix- 



ture of scalar and vector potentials of barrier type (13 1: 1) very strong (at small 



a) dependence of F on the binding energy of tunneling fermion and on the mixing 
coefficient A; 2) nonanalytic dependence of shift and width of level on the "force" 
(7 of scalar and vector long-range interactions. 

In conclusion, note that the version of SMSE considered above could seem to be 
rather artificial and having no relation to real problems. Let us emphasize thereupon 



that the Dirac equation with the potentials (13) at 1/2 < A ^ 1 can serve as the 



etalon equation for the relativistic theory of quasistationary states with the scalar- 
vector variant of interactions. As is known, now there are all grounds to consider 
that such interactions exist between composite objects (quarks and gluons) of QCD. 
Also, in the nuclear reactions of the tunnel type the peculiar features of a scalar- 
vector variant of interactions, that should be taken into account when calculating 
the penetration probability of tunneling fragments through potential barriers, can 
be fully manifested. 

Examples of strong influence of a scalar field on spectra of resonance states 
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of strongly interacting particles are given in the recent analysis (see, for example 
leviev^^ of experimental data from the Crystal Barrel Collaboration on the in-flight 
proton-antiproton annihilation into mesons of the final state below 2400 MeV. In 
particular, in the papei'^Sl the following assumption has been made. The existence 
of light cr-meson can be caused by singular behavior of quark-antiquark interaction 
~ 1/g^ at the small transmitted momentums. In the co-ordinate space it corresponds 
to linear increasing of the potential v{r) at large distances. 



Appendix A. Some notations 

The quantities introduced in Sec. 2 have the form 



a — b A± — c 

a — c 



^(^4, 5R = (l-.)(x^-.), (A.l) 



Xi-b"' y " (6 - d) 

xHb~c) 3H(6-c) ix^ -iy)if + 3c) , ^kjc"^ + c'' f + eg + h + l/c) 

(A.2) 



;(l-^) 



{b-cY 



f + ic 



3 (6-c)H 

4 ^ 



(A.3) 



1 r3(6-c)H2 23fi(3c2 + 2c/ + .g) 



43? 



(6-c) 



+ {b-c) [{\ + x^)v~ix^ 



+ H(/ + 3c)}, (A.4) 



/ 

(6 - c) 6c 
2 



, N, = [{b-\+){\+-c)]-\ Ar6 = [(&-A_)(A_-c)]-i, (A.5) 
' * H = x' (3 - 2i.) -(- (z/ - 2) . (A.6) 



(A+ -c)(A_ -c) \ 2(l-2A)cr 
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